Searches for gravitational waves from known pulsars with S5 

LIGO data 



B. P. Abbott^s, R. Abbott^s, F. Acemese^^'^^ R. Adhikari^s, R Ajith^, B. Allen^'^s, 
G. Allen^\ M. Alshourbagy^°''^ R. S. Amin^^, S. B. Anderson^^, W. G. Anderson^^, 

F. Antonucci^^", S. Aoudia^^", M. A. Arain*^^, M. Araya^®, H. Armandula^^, P. Armor^^, 
K. G. Arun^s, Y. Aso^^, S. Aston^^, R Astone^i^, P. Aufmuth^^, C. Aulbert^, S. Babak\ 
P. Baker3^ G. Ballardin^\ S. Ballmer^s, C. Barker^^, D. Barker^^, F. Barone^*"^ B. Barr^' 
P. Barriga^"^, L. Barsotti^i, M. Barsuglia^ M. A. Barton^^, I. Bartos^^, R. Bassiri^'', 

M. Bastarrika*'^, Th. S. Bauer^°", B. Behnke\ M. Beker^°, M. Benacquista5^ 
J. Betzwieser^s, P. T. Beyersdorf^^ S. Bigotta^O'^^ 1. A. Bilenko^^ G. Billingsley^^ , 

S. Birindelli^^a^ Biswas^^ M. A. Bizouard^^, E. Black^^, J. K. Blackburn^s, 
L. Blackburn^i, D. Blair^^ B. Bland^^, C. Boccaral^ T. P. Bodiya^i, L. Bogue^o, 
F. Bondu^26^ L Bonelli20«^ R. Bork^^, V. Boschi^s, S. Bose^^ L. Bosil9^ S. Braccini^o^ 

C. Bradaschia^°°, P. R. Brady^^ V. B. Braginsky^^^ J. E. Brau^^ D. O. Bridges^", 
A. Brillef^^", M. Brinkmann^, V. Brisson^^, C. Van Den Broeck*, A. F. Brooks^*, 

D. A. Brown^^ A. Brummit^^ G. Brunet3\ R. Budzynski^^^ T. Bulik^^'^'^, A. BuUington^^ 

H. J. Bulten^°'*^ A. Buonanno6^ O. Burmeister^, D. Buskulic^^, R. L. Byer^^ 
L. Cadonati6^ G. Cagnoli^6^ E. CallonilS'^^ J. B. Camp^^, E. Campagna^^'^^ 
J. Cannizzo^®, K. C. Cannon^^, B. Canuel^\ J. Cao^'^, F. Carbognani^^, L. Cardenas^^, 
S. Caride^^ G. Castaldi^\ S. Caudill^^, M. Cavaglia5^ F. Cavalier^s, R. Cavalierii\ 
G. Cella^o^ C. Cepeda^^, E. Cesarinii6'=, T. Chalermsongsak^s, E. Chalkley6^ 
P. Charlton^^ E. Chassande-Mottin^ S. Chatterji^^, S. Chelkowski^^^ Y. Chen^'^ 

A. Chincarinil^ N. Christensen^ C. T. Y. Chung5^ D. Clark5\ J. Clark^ J. H. Clayton^^ 
F. Cleva^2a^ ^ Coccia22'^^ T. Cokelaer^ C. N. Colacino^3,2o^ j_ Colas", A. Co\W^''\ 

M. Colombini^i^ R. Conte^^^ D. Cook^^, T. R. C. Corbitt^i, C. Corda20'^^ N. Cornish^e, 
A. Corsi^i'^^ J.-P. Coulon^2a^ Coward^^ D. C. Coyness, J. D. E. Creighton^^ 
T. D. Creighton5^ A. M. Cruise'^^^ R. M. Culter^^ A. Cumming6^ L. Cunningham^^, 

E. Cuoco", S. L. Danilishin^^, S. D'Antonio^^'^, K. Danzmann2'27^ A. Dari^9'^^ V. Dattilo^- 

B. Daudert^s, M. Davier^s, G. Davies^ E. J. Daw5^ R. Dayi\ R. De Rosa^S''^ D. DeBra^^ 
J. Degallaix^, M. del Prete^O'^'^, V. Dergachev^^ S. Desai^^^ R. DeSalvo^^, S. Dhurandhar^^, 

L. Di Fiore^s^ A. Di Ueto'^°''\ M. Di Paolo Emilio^^'^'^, A. Di Virgilio^O", M. Dfaz^^ 

A. Dietz8'26^ F. Donovan^i, K. L. Dooley^^^ E. E. Doomes^°, M. Drago^^'^'^, 
R. W. P. Drever^ J. Dueck^, I. Duke^S J.-C. Dumas^^ J. G. Dwyer^^, C. Echols2^ 
M. Edgar64, A. Effler^^, P. Elirens^^, E. Espinoza^^, T. Etzel^^, M. Evans^^, T. Evans^^, V. 
Fafone^^aft^ g Fairhurst^ Y. Faltas^^ Y. Fan^^, D. Fazios, H. Fehrmann^, I. Ferrante^"'^'', F 
Fidecaro20'^^ L. S. Finn^^ I. Fiori^^ R. Flaminio^^ K. Flasch^^ S. Foley3\ C. Forrest^°, 
N. Fotopoulos^^, J.-D. Fournier^^", J. Franc^^, A. Franzen^'', S. Frasca^^"^, F. Frasconi^"'', 



-2- 



M. Frede^, M. Frei^^ Z. Frei^^ A. Freise^^^ R. Frey'^^ T. Fricke^^, P. FritscheP\ 
V. V. Frolov^o, M. Fyffe^o, V. Galdi^\ L. Gammaitonil9'^^ J. A. Garofoli^^^ F. Garufil8"^ 

G. Gemme^'', E. Genin^\ A. Gennai^'^'', I. Gholami^, J. A. Giaime^^'^", S. Giampanis^, 
K. D. Giardina^o, A. Giazotto^o^ K. Goda^i, E. Goetz^^, L. M. Goggin^^^ G. Gonzalezes, 
M. L. Gorodetsky37, S. GoBler^-^o, R. Gouaty^^, M. Granata^ V. Granata^^, A. Grant6^ 

S. Gras^^ C. Gray^^, M. Gray^ R. J. S. Greenhalgh4^ A. M. Gretarsson^^, C. Greverie^^a^ 
F. Grimaldi^i, R. Grosso5^ H. Grote^, S. Grunewald\ M. Guenther^^, G. Guidi^^'^^ 
E. K. Gustafson^s, R. Gustafson^^, B. Hage^^, J. M. Hallam^^ D. Hammer^^^ 
G. D. Hammond^^, C. Hanna^^, J. Hanson^", J. Harms^^, G. M. Harry^^ I. W. Harry^, 

E. D. Harstad6^ K. Haughian6^ K. Hayama5^ J. Heefner^^, H. Heitmann^^ P. Hello^^, 

I. S. Heng64, A. HeptonstalF, M. Hewitson^, S. Hild'^^^ E. Hirose^^^ D. Hoak^o, 
K. A. Hodge^s, K. Holt^o, D. J. Hosken^i, J. Hough^^, D. Hoyland^^ D. Hueti\ 
B. Hughey3\ S. H. Huttner6^ D. R. Ingram^^, T. Isogai^ M. Ito'^^ A. Ivanov^^, 
P. Jaranowski^^*^, B. Johnson^^, W. W. Johnson'^'^, D. I. Jones^^, G. Jones^, R. Jones^^, 
L. Sancho de la Jordana^^, L. Ju^^, P. Kalmus^^, V. Kalogera^^, S. Kandhasamy^^, 
J. Kanner^^, D. Kasprzyk^^, E. Katsavounidis^^, K. Kawabe^^, S. Kawamura^^, 

F. Kawazoe^, W. Kells^^, D. G. KeppeP, A. Khalaidovski^, F. Y. Khalili^^^ R. Khan^o, 
E. Khazanov^e, P. King^^, J. S. KisseP^, S. Klimenko^^, K. Kokeyama^^, V. Kondrashov^^, 

R. Kopparapu^e, S. Koranda''^, I. Kowalska'^^'^, D. Kozak^^, B. Krislinan^, A. Krolak''^"-'^, 

R. Kumar''^ P. Kwee^^ P. La Pennai\ P. K. Lam^ M. Landry^^, B. Lantz5\ 
A. Lazzarini^*^, H. Lei^^, M. Lei^^, N. Leindecker^\ I. Leonor^^, N. Leroy^^, N. Letendre^^, 
C. Li^ H. R E. Lindquist^^, T. B. Littenberg^e, N. A. Lockerbie^^ D. Lodhia^^ 

M. Longo^\ M. LorenziniiS", V. Loriette^^, M. Lormand^^, G. Losurdo^'^'', P. Lu^\ 
M. Lubinski^s, A. LucianettF, H. Liick^'^^, B. Machenschalk\ M. Maclnnis^^ J.-M. 
Mackowski32, M. Ma geswaran , K. Mailand^^ E. Majorana^i", N. Man^^", I. Mandel'^^ 
V. Mandic6^ M. Mantovani^o^, F. Marchesoni^^'^, F. Marion^^, S. Marka^°, Z. Marka^^, 
A. Markosyan^\ J. Markowitz^i, E. Maros2^ J. Marque", F. Martelli^*'"^ I. W. Martin^^, 
R. M. Martin^^ J. N. Marx^^ K. Mason3\ A. Masserot^^ F. Maticharde^ L. Matone^°, 
R. A. Matzner^^, N. Mavalvala^i, R. McCarthy's, D. E. McClelland^ S. C. McGuire^^, 
M. McHugh3^ G. McIntyre2^ D. J. A. McKechan^ K. McKenzie^ M. Mehmet^ 
A. Melatos^^ A. C. Melissinos^°, G. MendelP^, D. F. Menendez^^^ F. Menzinger", 
R. A. Mercer^^ S. Meshkov'^, C. Messenger', M. S. Meyer^o, C. MicheF, L. Milano^^"^ 
J. Miller^^ J. Minelli^^ Y. Minenkov''"', Y. Mino^ V. P. Mitrofanov^^ G. Mitselmakher^^ 
R. Mittleman^i, O. Miyakawa'^ B. Moe^^ M. Mohan", S. D. Mohanty^^ 
S. R. P. Mohapatra^^ J. Moreau^^ G. Moreno'^ N. Morgado^', A. Morgia''''^ 
T. Morioka3^ K. Mors', S. Mosca^S''^ V. Moscatelli'^'^, K. Mossavi', B. Mours'^ 
C. MowLowry^, G. Mueller'^^, D. Muhammad^", H. zur Miihlen'^, S. Mukherjee^^, 

H. Mukhopadhyay'^ A. Mullavey^ H. Miiller-Ebhardt', J. Munch6\ P. G. Murray^^ 



- 3 - 



E. Myerses, J. Myers^^, T. Nash^^, J. Nelson6^ I. Neri^^ab^ q Newton6^ A. Nishizawa^^, 
F. Nocerai\ K. Numata^^, E. Ochsner6^ J. 0'Dell^^ G. H. Ogin^s, B. O'Reilly^o, 
R. 0'Shaughnessy^^ D. J. Ottaway^\ R. S. Ottens6^ H. Overmier^o, B. J. Owen^^^ 

G. Pagliaroli22'^'^, C. Palomba^i^ Y. Pan^s, C. Pankow^^, F. Paoletti20'^'i\ M. A. Papal'^^ 
V. Parameshwaraiah^^ , S. Pardi^^"*, A. Pasqualetti^^, R. Passaquieti^""^^, D. Passuello^°", 

P. PateP, M. Pedraza^s, S. Penn^^ A. Perreca^^^ G. PersichettilS''^ M. Pichot^^a^ 
F. Piergiovanni^^'^^ V. Pierro^\ M. Pietka^^^ L. Pinard^^ I. M. Pinto^\ M. Pitkin^^, 

H. J. Pletsch^, M. V. Plissi6^ R. Poggiani20'^^ F. Postiglione^^^ M. Prato^^ M. Principe^^ 
R. Prix^, G. A. ProdP"^ L. Prokhorov^^ O. Puncken^, M. Punturo^^^ P. Puppo2^^ 

V. Quetschke^^^ F. J. Raab^^, O. Rabaste^ D. S. Rabeling4°"^ H. Radkins^^, P. Raffai^^ 
Z. Raics^°, N. Rainer^, M. Rakhmanov5^ P. Rapagnani^^'^'', V. Raymond^\ V. Re^3"^ 
C. M. Reed29, T. Reed^^, T. Regimbau^^a^ ^ Rehbein^, S. Reid6^ D. H. Reitze^^^ 
F. Ricci^i'^^ R. Riesen^o, K. Riles^^ B. Rivera^^, P. Roberts^ N. A. Robertson^s-^^, 
F. Robinet^s, C. Robinson^ E. L. Robinson\ A. Rocchi^^", S. Roddy^", L. Rolland^^, 
J. Rollins^", J. D. Romano^^ R. Romano^^'^^ J. H. Romie^^, D. Rosinska^^^''^, C. Rover^, 
S. Rowan64, Rudiger^, P. Ruggi^\ P. Russell^^, K. Ryan^^, S. Sakata^^, F. Salemi4=^"^ 
V. Sandberg^^, V. Sannibale^^, L. Santamaria^, S. Saraf^^, P. Sarin'^^, B. Sassolas^^, 
B. S. Sathyaprakash^ S. Sato^^ M. Satterthwaite^ P. R. Saulson^^ R. Savage^^, P. Savov'', 
M. Scanlan^^, R. Schilling^, R. Schnabel^, R. Schofield6^ B. Schulz^, B. F. Schutz^'^ 
P. Schwinberg29, J. Scott^^, S. M. Scott^, A. C. Searle^^, B. Sears^^, F. Seifert^, 
D. Sellers^°, A. S. Sengupta^*^, D. Sentenac^\ A. Sergeev^^, B. Shapiro^^ P. Shawhan^^, 

D. H. Shoemaker^i, A. Sibley^o, X. Siemens^^ D. Sigg^^, S. Smha^\ A. M. Sintes'^o, 
B. J. J. Slagmolen^ J. Slutsky^^, M. V. van der Sluys^\ J. R. Smith^^^ M. R. Smithes, 
N. D. Smith^i, K. Somiya^ B. Sorazu6^ A. Stein^i, L. C. Stein^i, S. Steplewski^^ 
A. Stochino^s, R. Stone^^ K. A. Strain6^ S. Strigin^^, A. Stroeer^s, R. Sturanii^'^^ 
A. L. Stuver^o, T. Z. Summerscales^, K. -X. Sun^\ M. Sung^^, P. J. Sutton^ B. Swinkels", 
G. P. Szokolyi^^ D. Talukder^^ L. Tang5^ D. B. Tanner^^, S. P. Tarabrin^^, J. R. Taylor^, 

R. Taylor^s, R. Terenzi^^'^^, J. Thacker=^°, K. A. Thorne^", K. S. Thorne^ A. Thiiring^^, 
K. V. Tokmakov^^, A. Toncelli^O'^'', M. Tonelli^O"^ C. Torres=^°, C. Torrie^^, E. Tournefier^e, 

F. Travassol*^"^ G. Traylor^o, M. Trias''°, J. Trummer^e, D. Ugolini^^ J. mmen^\ 
K. Urbanek^i, H. Vahlbruch^^, G. Vajente^O"'', M. Vallisneri^ J. F. J. van den Brand^°'^^ S. 
van der Putten^°", S. Vass^^, R. Vaulin^^ M. Vavoulidis^^ , A. Vecchio^^^ G. Vedovato^^c^ 
A. A. van VeggeP, J. Veitch^^^ P. Veitch^S C. Veltkamp^, D. Verkindt^^, F. Vetrano^''"^ 
A. Vicere^S'^^ A. Villains, J.-Y. Vinet^^a^ ^ Vocca^^", C. Vorvick^^, S. P. Vyachanin^^, 
S. J. Waldman^i, L. Wallace^^, R. L. Ward^^, M. Was^^, A. Weidner^, M. Weinert^, 
A. J. Weinstein^s, R. Weiss^^, L. Wen^■^^ S. Wen^^, K. Wette^ J. T. Whelan^•4^ 
S. E. Whitcomb2^ B. F. Whiting^^^ C. Wilkinson^^, P. A. Willems^s, H. R. Williams^^^ 
L. Wmiams*53^ B. Willke^'^^, I. Wilmut^^ L. Winkelmann^ , W. Winkler^, C. C. Wipf^^, 



-4- 



A. G. Wiseman^^ G. Woan'5^ R. Wooley^o, J. Worden^^, W. Wu^^^ 1. Yakushin^o, 
H. Yamamoto^s, Z. Yan^^ S. Yoshida^^, M. Yvert^e, M. Zanolin^^ J. Zhang^^ L. Zhang^^, 
C. Zhao^^ N. Zotov^^ M. E. Zucker^i, J. Zweizig^^ 

The LIGO Scientific Collaboration & The Virgo Collaboration 



- 5 - 



^Albert-Einstein-Institut, Max-Plaiick-Institut fur Gravitationsphysik, D-14476 Golm, Germany 
^Albert-Einstein-Institut, Max-Planck-Institut fiir Gravitationsphysik, D-30167 Hannover, Germany 
■^Andrews University, Berrien Springs, MI 49104 USA 

^AstroParticulc ct Cosmologic (AFC), CNRS: UMR7164-IN2P3-Observatoirc dc Paris- Universitc Denis 
Diderot-Paris VII - CEA : DSM/IRFU 

■^Australian National University, Canberra, 0200, Australia 

^California Institute of Technology, Pasadena, CA 91125, USA 

^Caltech-CaRT, Pasadena, CA 91125, USA 

8 Cardiff University, Cardiff, CF24 3AA, United Kingdom 

^Carleton College, Northfield, MN 55057, USA 
''''Charles Sturt University, Wagga Wagga, NSW 2678, Australia 
lOColumbia University, New York, NY 10027, USA 

^^European Gravitational Observatory (EGO), 1-56021 Cascina (Pi), Italy 
i^Embry-Riddle Aeronautical University, Prescott, AZ 86301 USA 
^^Eotvos University, ELTE 1053 Budapest, Hungary 
i-^ESPCI, CNRS, F-75005 Paris, France 

i^Hobart and William Smith Colleges, Geneva, NY 14456, USA 

^^INFN, Sezione di Firenze, 1-50019 Sesto Fiorentino'^; Universita degli Studi di Firenze, 1-50121'', Firenze; 
Universita degli Studi di Urbino 'Carlo Bo', 1-61029 Urbino^, Italy 

^"^INFN, Sezione di Genova; 1-16146 Genova, Italy 

^^INFN, sezione di Napoli °; Universita di Napoli 'Federico H''' Complesso Universitario di Monte S.Angelo, 
1-80126 Napoh; Universita di Salerno, Fiseiano, 1-84084 Salerno'^, Italy 

^^INFN, Sezione di Perugia"; Universita di Perugia'', 1-6123 Perugia, Italy 

20INFN, Sezione di Pisa*^; Universita di Pisa''; 1-56127 Pisa; Universita di Siena, 1-53100 Siena^ Italy 
^^INFN, Sezione di Roma"; Universita 'La Sapienza''', 1-00185 Roma, Italy 

^^INFN, Sezione di Roma Tor Vergata"; Universita di Roma Tor Vergata'', Istituto di Fisica dello Spazio 
Interplanetario (IFSI) INAF^ 1-00133 Roma; Universita deU'Aquila, T67100 L'Aquila'', Italy 

^•^Institute of Applied Physics, Nizhny Novgorod, 603950, Russia 

^"'inter-University Centre for Astronomy and Astrophysics, Pune - 411007, India 

25LAL, Universite Paris-Sud, IN2P3/CNRS, F-91898 Orsay, France 

^^Laboratoire d'Annecy-le-Vieux de Physique des Particules (LAPP), IN2P3/CNRS, Universite de Savoie, 



- 6 - 



F- 74941 Annecy-le-Vieux, France 

^^Leibniz Universitat Hannover, D-30167 Hannover, Germany 

28LIGO - California Institute of Technology, Pasadena, CA 91125, USA 

29LIGO - Hanford Observatory, Richland, WA 99352, USA 

3°LIG0 - Livingston Observatory, Livingston, LA 70754, USA 

31LIG0 - Massachusetts Institute of Technology, Cambridge, MA 02139, USA 

^^Laboratoire des Materiaux Avances (LMA), IN2P3/CNRS, F-69622 Villeurbanne, Lyon, France 

■^■^Louisiana State University, Baton Rouge, LA 70803, USA 

^^Louisiana Tech University, Ruston, LA 71272, USA 

^^Loyola University, New Orleans, LA 70118, USA 

^•^Montana State University, Bozeman, MT 59717, USA 

•^^Moscow State University, Moscow, 119992, Russia 

38NASA/Goddard Space Flight Center, Greenbelt, MD 20771, USA 

■^^National Astronomical Observatory of Japan, Tokyo 181-8588, Japan 

"ONikhef, National Institute for Subatomic Physics, P.O. Box 41882, 1009 DB Amsterdam, The 
Netherlands"; VU University Amsterdam, De Boelelaan 1081, 1081 HV Amsterdam, The Netherlands'' 

^^Northwestern University, Evanston, IL 60208, USA 

^^Departement Artemis, Observatoire de la Cote d'Azur, CNRS, F-06304 Nice Institut de Physique de 
Rennes, CNRS, Universite de Rennes 1, 35042 Rennes ^; France 

"^^INFN, Gruppo CoUegato di Trento'' and Universita di Trento^ L38050 Povo, Trento, Italy; INFN, 
Sezione di Padova^ and Universita di Padova'', 1-35131 Padova, Italy 

44IM-PAN 00-956 Warsaw"^; Warsaw Univ. 00-681^ Astro. Obs. Warsaw Univ. 00-478^=; CAMK-PAM 
00-716 Warsaw"*; Bialystok Univ. 15-424^ IPJ 05-400 Swierk-Otwock^; Inst, of Astronomy 65-265 Zielona 
Gora ^ , Poland 

"^Rochester Institute of Technology, Rochester, NY 14623, USA 

^^Rutherford Appleton Laboratory, HSIC, Chilton, Didcot, Oxon 0X11 OQX United Kingdom 
■*^San Jose State University, San Jose, CA 95192, USA 

Sonoma State University, Rohnert Park, CA 94928, USA 
"^^ Southeastern Louisiana University, Hammond, LA 70402, USA 
soSouthcrn University and A&M College, Baton Rouge, LA 70813, USA 
siStanford University, Stanford, CA 94305, USA 
^^Syracuse University, Syracuse, NY 13244, USA 



- 7- 



S. Begm86'S^ A. Corongiu^^ N. D'Amico^^.s!^ p_ q preire^^'^^, J. W. T. Hessels^^'^^ 
G. B. Hobbs^o, M. KramerS^ A. G. Lyne^^ R. N. Manchester^^, F. E. MarshallS^ 
J. Middleditch^^ A. PossentP^ S. M. Ransom^^, I. H. Stairs^^, and B. Stappers^^ 



^^Thc Pennsylvania State University, University Park, PA 16802, USA 
5'*The University of Melbourne, Parkville VIC 3010, Australia 
'"'^The University of Mississippi, University, MS 38677, USA 
56The University of Sheffield, Sheffield SIO 2TN, United Kingdom 
■"The University of Texas at Austin, Austin, TX 78712, USA 

'"'^The University of Texas at Brownsville and Texas Southmost College, Brownsville, TX 78520, USA 

s^Trinity University, San Antonio, TX 78212, USA 

^^Universitat de les Illes Balears, E-07122 Palma de Mallorca, Spain 

^^University of Adelaide, Adelaide, SA 5005, Australia 

^^University of Birmingham, Birmingham, B15 2TT, United Kingdom 

s^University of Florida, Gainesville, FL 32611, USA 

^^University of Glasgow, Glasgow, G12 8QQ, United Kingdom 

e^University of Maryland, College Park, MD 20742 USA 

'^''Umversity of Massachusetts - Amherst, Amherst, MA 01003, USA 

s^University of Michigan, Ann Arbor, MI 48109, USA 

•^^University of Minnesota, Minneapohs, MN 55455, USA 

e^University of Oregon, Eugene, OR 97403, USA 

^°University of Rochester, Rochester, NY 14627, USA 

^^University of Sannio at Bcnevento, 1-82100 Benevento, Italy 

''^University of Southampton, Southampton, S017 IBJ, United Kingdom 

^■^University of Strathclyde, Glasgow, Gl IXQ, United Kingdom 

^^Univcrsity of Western Australia, Crawley, WA 6009, Australia 

^^University of Wisconsin-Milwaukee, Milwaukee, WI 53201, USA 

^'^Washington State University, Pullman, WA 99164, USA 



- 8 - 



ABSTRACT 



We present a search for gravitational waves from 116 known millisecond and 
young pulsars using data from the fifth science run of the LIGO detectors. For 
this search ephemerides overlapping the run period were obtained for all pulsars 
using radio and X-ray observations. We demonstrate an updated search method 
that allows for small uncertainties in the pulsar phase parameters to be included 
in the search. We report no signal detection from any of the targets and therefore 
interpret our results as upper limits on the gravitational wave signal strength. 
The most interesting limits are those for young pulsars. We present updated 
limits on gravitational radiation from the Crab pulsar, where the measured limit 
is now a factor of seven below the spin-down limit. This limits the power radiated 
via gravitational waves to be less than ~2% of the available spin-down power. For 
the X-ray pulsar J0537— 6910 we reach the spin-down limit under the assumption 
that any gravitational wave signal from it stays phase locked to the X-ray pulses 
over timing glitches, and for pulsars J1913-I-1011 and J1952+3252 we are only 
a factor of a few above the spin-down limit. Of the recycled millisecond pulsars 
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several of the measured upper limits are only about an order of magnitude above 
their spin-down limits. For these our best (lowest) upper limit on gravitational 
wave amplitude is 2.3x10"^^ for J1603— 7202 and our best (lowest) limit on the 
inferred pulsar ellipticity is 7.0 x 10"^ for J2124-3358. 

Subject headings: gravitational waves - pulsars: general 



1. Introduction 

Within our Galaxy some of the best targets for gravitational wave searches in the sen- 
sitive frequency band of current interferometric gravitational wave detectors (~40-2000 Hz) 
are millisecond and young pulsars. There are currently just over 200 known pulsars with 
spin frequencies greater than 20 Hz, which therefore are within this band. In this paper we 
describe the latest results from the ongoing search for gravitational waves from these known 
pulsars using data from the Laser Interferometric Gravitational- Wave Observatory (LIGO). 
As this search looks for objects with known positions and spin-evolutions it can use long 
time spans of data in a fully coherent way to dig deeply into the detector noise. Here we use 
data from the entire two-year run of the three LIGO detecto rs, entitled Science R un 5 (S5), 



during which the detectors reached their design sensitivities ( lAbbott et al.ll2009bl ). This run 
started on 2005 November 4 and ended on 2007 October 1. The detectors (the 4 km and 2 km 
detectors at LIGO Hanford Observatory, HI and H2, and the 4 km detector at the LIGO Liv- 
ingston Observatory, LI) had duty factor s of 78% for HI, 79% for H2, and 66% for LI. Th e 



GEO 600 detector also participated in S5 (iGrote fc the LIGO Scientific Collaborationll2008l ) 



but at lower sensitivities that meant it was not able to enhance this search. T he Virgo detec- 



tor a lso had data overlapping with S5 during Virgo Science Run 1 (VSRl) ( lAcernese et al. 



20081 ). However this was also generally at a lower sensitivity than the LIGO detectors and 
had an observation time of only about 4 months, meaning that no significant sensitivity 
improvements could be made by including this data. Due to its multi-stage seismic isola- 
tion system Virgo does have better sensitivity than the LIGO detectors below about 40 Hz, 
opening the possibility of searching for more young pulsars, including the Vela pulsar. These 
lower frequency searches will be explored more in the future. 

This search assumes that the pulsars are triaxial stars emitting gravitational waves at 
precisely twice their observed spin frequencies, i.e. the emission mechanism is an £ = m = 2 
quadrupole, and that gravitational waves are phase-locked with the electromagnetic signal. 
We use the so-called spin-down limit on strain tensor amplitude h^f^ as a sensitivity target 
for each pulsar in our analysis. This can be calculated, by assuming that the observed spin- 
down rate of a pulsar is entirely due to energy loss through gravitational radiation from an 



-lo- 



rn 



2 quadrupole, as 



8.06xl0-i%8r-pi^(|z>|/z/)i/2, 



(1) 



where J38 is the pulsar's principal moment of inertia (Izz) in units of lO'^^kgm^, r^pc is 
the pulsar distance in kpc, u is the spin-frequency in Hz, and i> is the spin-down rate in 
Hzs~^. Due to uncertainties in Izz and r, h^^ is typically uncertain by about a factor 2. 
Part of this is due to the uncertainty in Izz which, though predicted to lie roughly in the 
range 1-3 x 10^^ kg m^, has not been measured for any neutron star; and the best (though 
still uncertain) pr ospect is star A of th e double pulsar system J0737-3039 with 20 years' 
more observation (IKramer fc Wexll2009l ). Distance estimates based on dispersion measure 
can also be w rong by a factor 2-3, as confirmed by recent parallax observations of the 
double pulsar jDeller et aliboOQh . For pulsars with measured braking indices, n = vvjiP'^ 
the assumption that spin-down is dominated by gravitational wave emission is known to be 
false (the braking index for quadrupolar gravitational wave emission should be 5, but all 
measured n's are less than 3) and a stricter indirect limit on gravit ational wave ein ission can 
be set. A phenomenological investigation of some young pulsars (jPalombal l2000l ) indicates 
that this lim it is lo wer t han by a f actor 2.5 or more, depending on the pulsar. See 
Abbott et al.l (120071 ) and I Abbott et al.l (120081 ) for more discussion of the uncertainties in 
indirect limits. Recycled millisecond pulsars have intrinsically small spin-downs, so for the 
majority of pulsars in our search these spin-down limits will be well below our current 
sensitivities, making detection unlikely. However, our search also covers four young pulsars 
with large spin-down luminosities, and for these we can potentially beat or reach their spin- 
down limits using current data. 

The LIGO band covers the fastest (highest-z/) known pulsars, and the quadrupole for- 
mula for strain tensor amplitude 



/lo = 4.2 X lO^^^i^i^oo^as^-erkp, 



pc 



(2) 



indicates that these pulsars are the best gravitational wave emitters for a given equatorial 
ellipticity e = {I^x — Iyy)/Izz (here i^ioo = i^/(100 Hz) and £_6 = e/lO"^). The pulsars with 
high spin-downs are almost all less than ~ 10^ years old. Usually this is interpreted as greater 
electromagnetic activity (including particle winds) in younger objects, but it could also mean 
that they are more active in gravitational wave emission. This is plausible o n theoretica l 
grounds too. Strong internal magnetic fields may cause significant elhpticities (ICutlerl 120021) 
which would then decay as the field decays or otherwise changes (iGoldreich fc Reisenegger 
I992I ). The initial crust may be asymmetric if it forms on a time scale on which the neutron 
star is still perturbed by its violent fo rmation and aftermath, including a p ossible lengthy 
perturbation due to the fluid r-modes (iLindblom et al.ll2000l : IWu et al.ll200ll ). and asymme- 
tries may slowly relax due to mechanisms such as viscoelastic creep. Also the fluid r-modes 
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may remain unstable to gravitational wave emission for up to a few thousand years after 
the neutron star's birth, depending on its composition, viscosity, and initial spin frequency 
( lOwen et al.lll998l : iBondarescu et al.ll2009[ ). Such r-modes are expected to have a gravita- 
tional wave frequency about 4/3 the spin frequency. However, we do not report on r-mode 
searches in this paper. 



1.1. Previous analyses 



The first search for gravitational waves from a known pulsar using LIGO and GEO 600 
data came from the first science run (SI) in 2002 September. This targeted just one pulsar 
i n the approximatel y one weeks worth of data - the then fastest known pulsar J1939+2134 
( lAbbott et al.l 120041 ). Data from LIGO's second science run (S2), which spanned from 2003 
February to 2003 April, was used to search for 28 isolated pulsars (i.e. those not in binary 
systems) ( lAbbott et al.ll2005l ). The last search for gravitational waves from multiple known 
pulsars using LIGO data combined data from the third and fourt h science runs and had 78 
targets, including isolated pulsars and those in binary systems (lAbbott et al.l 120071 ). The 
best (lowest), 95% degree-of-belief, upper limit on gravitational wave amplitude obtained 
from the search was /iq^^" = 2.6 x 10~^^ for J1603— 7202, and the best (smallest) limit on 
ellipticity was just under 10~^ for J2124— 3358. The data run used in this paper is almost 
an order of magnitude longer, and has a best strain noise amplitude around a factor of two 
smaller, than that used in the best previous search. 

We have also previously sear ched the first nine months of S5 data for a signal from the 
Crab pulsar (lAbbott et al.ll2008l ). That analysis used two methods to search for a signal: 
one in which the signal was assumed to be precisely phase-locked with the electromagnetic 
signal, and another which searched a small range of frequencies and frequency derivatives 
around the electromagnetic parameters. The time span of data analysed was dictated by 
a timing glitch in the pulsar on 2006 August 23, which was used as the end point of the 
analysis. In that search the spin-down limit for the Crab pulsar was beaten for the first 
time (indeed it was the first time a spin-down limit had been reached for any pulsar), with 
a best limit of /ig^^° = 2.7x10"^^, or slightly below one-fifth of the spin-down limit. This 
allowed the total power radiated in gravitational waves to be constrained to less than 4% 
of the spin- down power. We ha ve since discovered an error in the signal template used for 
the search (lAbbott et al.ll2009al ). We have re-analysed the data and find a new upper limit 
based on the early S5 data alone at the higher value shown in Table [3l along with the smaller 
upper limit based on the full S5 data. 



For this analysis we have approximately 525 days of HI data, 532 days of H2 data and 



-VI- 



437 days of LI data. This is using all data flagged as science mode during the run (i.e. taken 
when the detector is locked in its operating condition on the dark fringe of the interference 
pattern, and relatively stable), except data one minute prior to loss of lock, during which 
time it is often seen to become more noisy. 



1.2. Electromagnetic observations 



The radio pulsar parameters used for our searches are based on ongoing radio pulsar 
monitoring programs, using data from the Jodrell Bank Observatory (JBO), the NRAO 
100 m Green Bank Telescope (GBT) and the Parkes radio telescope of the Australia Telescope 
National Facility. We used radio data coincident with the S5 run as these would reliably 
represent the pulsars' actual phase evolution during ou r searches. We obtained data for 44 
pulsars from JBO (including the Crab pulsar ephemeris, iLyne et al.l (Il993l . l2009l )). 39 pulsars 
within the Terzan 5 and M28 globular clusters from G BT, and 47 from P arkes, including 
pulsars timed as part of the Parkes Pulsar Timing Array (jManchesterll2008l ). For 15 of these 
pulsars there were observations from more than one site, making a total of 115 radio pulsars 
in the analysis (see Table [1] for list of the pulsars, including the observatory and time span of 
the observations). For the pulsars observed at JBO and Parkes we have obtained parameters 
fit to data overlapping with the entire S5 run. For the majority of pulsars observed at GBT 
the parameters have been fit to data overlapping approximately the first quarter of S5. 

Pulsars generally exhibit timing noise on long time scales. Over tens of years this can 
cause correlations in the pulse time of arrivals which can give systematic errors in the pa- 
rameter fits produced, by the standard pulsar timing packag e TEMPoF^. of orde r 2-10 times 
the regular errors that TEMPO assigns to each parameter (jVerbiest et al.ll2008l ). depending 
on the amplitude of the noise. For our pulsars, with relatively short observation periods 
of around two years, the long-term timing noise variations should be largely folded in to 
the parameter fitting, leaving approximately white uncorrelated residuals. Also millisecond 
pulsars, in general, have intrinsically low levels of timing noise, showing comparatively white 
residuals. This should mean that the errors produced by TEMPO are approximately the 
true \a errors on the fitted values. 



The regular pulse timing observations of the Crab pulsar (ILyne et al.l Il993l . |2009| ) in- 
dicate that the 2006 August 23 glitch was the only glitch during the S5 run. One other 
radio pulsar, J1952-I-3252, was observed to glitch during the run (see §5.1.31 ) Independent 
ephemerides are available before and after each glitch. 



http : //www. atnf . csiro . au/research/pulsar/tempo/ 
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We include one pulsar in our analysis that is not observed as a radio pulsar. This is 
PSR J0537— 6910 in the Large Magellanic Cloud, for which only X-ray timings currently exist. 
Data for this source corn e from dedicated time on the Rossi X-ray Timing Explorer (RXTE) 
( iMiddleditch et al.l 120061 ) . giving ephemerides covering the whole of S5. These ephemerides 
comprise seven inter-glitch segments, each of which produces phase-stable timing solutions. 
The segments are separated by times when the pulsar was observed to glitch. Due to the 
complexity of the pulsar behaviour near glitches, which is not reflected in the simple model 
used to predict pulse times of arrival, sometimes up to ~ 30 days around them are not 
covered by the ephemerides. 



Gravitational wave search method 



The details of the search method are discussed in lDupuis &: Woaru (120051 ) and I Abbott et al 



(120071 ). but we will briefly review them here. Data from the gravitational wave detectors are 
heterodyned using twice the known electromagnetic phase evolution of each pulsar, which 
removes this rapidly varying component of the signal, leaving only the daily varying ampli- 
tude modulation caused by each detector's antenna response. Once heterodyned the (now 
complex) data are low-pass filtered at 0.25 Hz, and then heavily down-sampled, by averaging, 
from the original sample rate of 16 384 Hz to 1/60 Hz. Using these down-sampled data {B^, 
where k represents the k^^ sample) we perform parameter estimation over the signal model 
yk{ci) given the unknow n signal parameter s a. This is done by calculating the posterior 
probability distribution ( lAbbott et al.l 120071 1 



A4 



p{a\{B,}) cx n EW^'^-) - "^{yM}? + mSk} - "^{Vkia)})' x p{a), (3) 



where the first term on the right hand side is the likelihood (marginalised over the data 
variance, giving a Student 's-t-like distribution), p{a) is the prior distribution for a, M is 
the number of data segments into which the BkS have been cut (we assume stationarity of 
the data during each segment), rrij is the number of data points in the jth segment (with a 
maximum value of 30, i.e. we only assume stationarity for periods less than, or equal to, 30 
minutes in length), and n = Yl'j=i'^j- The assump tion of Gaussia nity and stationarity of 
the segments holds well for this analysis (see §4.5 of iDupuid (120041 ) for examples of cind 
Kolmogorov-Smirnov tests performed to assess these in previous analyses). 



We have previously (lAbbott et al.l |2004J . |2005| . 120071 ) performed parameter estimation 
over the four unknown gravitational wave signal parameters of amplitude ho, initial phase 0o, 
cosine of the orientation angle cost, and polarisation angle ip: giving a = {ho,(j)Q,cosL,ip}. 
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Priors on each parameter are set to be uniform over their allowed ranges, with the upper end 
of the range for Hq set empirically from the noise l evel of the data. We choose a uniforni prior 
on hn for co nsistency with our previous analyses (IDupuis fc Woanll2005l : lAbbott et al.ll2004l . 



2005 



20071 ) and to facilitate straightforward comparison of sensitivity. Extensive trials with 



software injections have shown this to be a very reasonable choice, returning a conservative 
(i.e. high) upper limit consistent with the data and any possible signal. 

Using a uniformly spaced grid on this four- dimensional parameter space the posterior is 
calculated at each point. To obtain a posterior for each individual parameter we marginalise 
over the three others. Using the marginalised posterior on ho we can set an upper limit 
by calculating the value that, integrating up from zero, bounds the required cumulative 
probability (which we have taken as 95%). We also combine the data from multiple detectors 
to give a joint posterior. To do this we simply take the product of the likelihoods for each 
detector and multiply this joint likelihood by the prior. This is possible due to the phase 
coherence between detectors. Again we can marginalise to produce posteriors for individual 
parameters. 

Below, in §2.1^ we discuss exploring and expanding this parameter space to more di- 
mensions using a Markov chain Monte Carlo (MCMC) technique. 



2.1. MCMC parameter search 

When high resolutions are needed it can be computationally time consuming to calculate 
the posterior over an entire grid as described above, and redundant areas of parameter space 
with very little probability are explored for a disproportionately large amount of time. A 
more efficient way to carry out such a search is with a Markov chain Monte Carlo (MCMC) 
technique, in which the parameter space is explored more efficiently and without spending 
much time in the areas with very low probability densities. 

An MCMC integration explores the parameter space by stepping from one position 
in parameter space to another, comparing the posterior probability of the two points and 
using a simple algorithm to determine whether the step should be accepted. If accepted 
it moves to that new position and repeats; if it is rejected it stays at the current position 
and repeats. Each iteration of the chain, whether it stays in the same position or not, is 
recorded and the amount of time the chain spends in a particular part of parameter space is 
directly proportional to the posterior probability density there. The new points are drawn 
randomly from a specific proposal distribution, often given by a multivariate Gaussian with 
a mean set as the current position, and a predefined covariance. For an efficient MCMC the 
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proposal distribution should reflect the underlying posterior it is sampling, but any proposal 
(that does not explicitly exclude the posterior), given enough time, will sample the posterior 
and deliver an accurate result. We use the Metropolis-Hastings (MH) algorithm to set the 
acceptance/rejection ratio. Given a current position MH accepts the new position Oj+i 
with probability 



where p{a\d) is the posterior value at a given data d, and q{ci\b) is the proposal distribution 
defining how we choose position a given a current position b. In our case we have symmet- 
ric proposal distributions, so g(aj+i|a)/g(aj|aj+i) = 1 and therefore only the ratio of the 
posteriors is needed. 

A well-tuned MCMC will efficiently explore the parameter space and generate chains 
that, in histogram form, give the marginalised posterior distribution for each parameter. 
Defining a good set of proposal distributions for the parameters in a has been done experi- 
mentally assuming that they are uncorrelated and therefore have independent distributions. 
(There are in fact correlations between the ho and cos l parameters and the 0o and ip pa- 
rameters, but in our studies these do not significantly alter the efficiency from assuming 
independent proposals.) The posterior distributions of these parameters will also generally 
not be Gaussian, especially in low SNR cases (which is the regime in which we expect to 
be), but a Gaussian proposal is easiest to implement and again does not appear to sig- 
nificantly affect the chain efficiency. We find that, for the angular parameters, Gaussian 
proposal distributions with standard deviations of an eighth the allowed parameter range 
(i.e. cr^p = 7r/4rad, (Jcost = 1/4 and = 7r/16rad) provide a good exploration of the pa- 
rameter space (as determined from the ratio of accepted to rejected jumps in the chain) for 
low SNR signals. We have performed many simulations comparing the output of the MCMC 
and grid-based searches, both on simulated noise and simulated signals, and both codes give 
results consistent to within a few percent. In these tests we find that the computational 
speed of the MCMC code is about three times faster than the grid-based code, although this 
can vary by tuning the codes. 

An MCMC integration may take time to converge on the bulk of the probability distri- 
bution to be sampled, especially if the chains start a long way in parameter space from the 
majority of the posterior probability. Chains are therefore allowed a burn-in phase, during 
which the positions in the chain are not recorded. For low SNR signals, where the signal 
amplitude is close to zero and the posteriors are reasonably broad, this burn-in time can be 
short. To aid the convergence we use simulated annealing in which a temperature parame- 
ter is used to fiatten the posterior during burn-in to help the chain explore the space more 
quickly. We do however use techniques to assess whether our chains have converged (see 
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Brooks fc Roberta (119981 ) for a good overview of convergence assessment tests for MCMCs.) 



We use two such tests: the Geweke test is used on individual chains to compare the means 
of two independent sections of the chain; and the Gelman and Rubins test is used to com- 
pare the variances between and within two or more separate chains. These tests are never 
absolute indicators of convergence, so each chain also has to be examined manually. The 
acceptance/rejection ratio of each chain is also looked at as another indicator of convergence. 
For all our results, discussed in §4] and ^ we have run these tests on the output chains, and 
conclude that all have converged. 



2.2. Adding phase parameters 

The heterodyne phase is calculated using the parameters measured from electromagnetic 
observations, which have associated errors. These errors could mean that the heterodyne 
phase will be offset, and drift away from, the true phase of the signal. Previously we have 
required that our data be heterodyned with a phase that was known to match the true 
signal phase over the course of the run to a few degrees. The criterion used to decide on 
whether to keep, or discard, a pulsar from the analysis was that there was no more than a 
30° drift of the elec tromagnetic phase from the signal phase over the course of a data run 



(lAbbott et all 120071 ) (i.e. if a signal was present t he phase drift would lead to a loss in SNR 



of less than about 15%.) In lAbbott et al.l (120071 ) this potential phase drift was calculated 



from the known uncertainties of the heterodyne phase parameters, but without taking into 
account the covariance between parameters, and as such was an over- conservative estimate, 
risking the possibility that some pulsars were excluded from the analysis unnecessarily. 

Rather than just setting an exclusion criterion for pulsars based on the potential phase 
mismatch (see Q we can instead search over the uncertainties in the phase parameters. This 
search can also be consistent with the, now provided, covariances on the phase parameters. 
For pulsars that have small mismatches over the run (that would have been included in the 
previous analysis), the extra search space allows these small uncertainties to be naturally 
folded (via marginalisation) into the final posteriors on the four main gravitational wave pa- 
rameters and our eventual upper limit. For pulsars with larger mismatches, which previously 
would have been excluded, this extra search space allows us to keep them in the analysis 
and again fold the phase parameter uncertainties into the final result. 

We can incorporate the potential phase error into the search by including the phase 
parameters, as an offset from the heterodyne phase parameters, in the parameter estimation 
and marginalising over them. A pulsar with, for example, associated errors on z/, z>, right 
ascension, declination, proper motion in both positional parameters and 5 orbital binary 
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parameters would add 11 extra parameters to the search. An MCMC is a practical way that 
allows us to search over these extra parameters. It also means that we make sure we cover 
enough of the parameter space so as not to miss a signal with a slight phase offset from the 
heterodyne values. Ex ampl es of MCMCs being used in a similar context can be found in 
Umstatter et al.l ( 12004| ) and IVeitch et al.l (120051 ) . However both these examples attempted 
to explore far greater parameter ranges than will be applied here. We do not attempt to 
use the MCMC as a parameter estimation tool for these extra parameters (i.e. due to our 
expected low SNR regime we would not expect to improve on the given uncertainties of the 
parameters), but just as a way of making sure we fully cover the desired parameter space, 
and fold the uncertainties into our final result without excessive computational cost. 

The number of extra parameters included in the search depends on how many parameters 
were varied to fit the ra dio data. When fit ting parameters with the pulsar timing packages 
TEMPO, or TEMP02 Jflobbs et al.lboOGh . certain values can be held fixed and others left 
free to vary, so errors will only be given on those allowed to vary. The uncertainties on those 
parameters that are fit will contain all the overall phase uncertainty. 



2.3. Setting up the MCMC 

In our search we have information on the likely range over which to explore each param- 
eter as given by the error from the fit (which we take as being the la value of a Gaussian 
around the best fit value), and the associated parameter covariance matrix. We use this 
information to set a prior on these parameters b given by a multivariate Gaussian 

pib) oc exp |-i(6 - boVC-\b - 6o)|, (5) 

with a covariance matrix C. For the vast majority of pulsars we expect that the uncertainties 
on the parameters are narrow enough that within their ranges they all give essentially the 
same phase model (see discussion of phase mismatch in §3])- In this case the posterior on 
the parameters should be dominated by this prior, therefore a good proposal distribution to 
efficiently explore the space with is the same multivariate Gaussian. 

An example of the posteriors produced when searching over additional parameters (in 
this case changes in declination and right ascension) can be seen in Figure [T] The figure 
shows the multivariate Gaussian used as priors on the two parameters and how the posterior 
is essentially identical to the prior (i.e. the data, which contains no signal, is adding no new 
information on those parameters, but the full prior space is being explored). We have assessed 
this technique on many simulations of noise and found that, as expected, the posteriors on the 
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additional phase parameters match the priors. We have also tested this technique on many 
software simulated signals (using real pulsar phase parameters) by offsetting the injected 
signal parameters from the the "best fit" values by amounts consistent with the parameter 
covariance matrix. In these tests we have been able to extract the parameters as expected. 




Fig. 1. — The posteriors and priors on offsets in declination 6 and right ascension a computed 
from an MCMC search for PSR J0407+1607 using a day of simulated data containing no 
signal. The covariance contour plot of the MCMC chains for the two parameters is shown 
and has a correlation coefficient of —0.93, which is identical to that of the multivariate 
Gaussian prior distribution used in this study. 
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2.4. Hardware injections 



During all LIGO science runs, except the first, fake pulsar signals have been injected 
into the detectors by direct actuation of the end test masses. These have provided end-to- 
end validation of the an alysis codes with coh erence between the different detectors. During 
S5, as with S3 and S4 ( lAbbott et al.l 120071 ). ten signals with different source parameters 
were injected into the detectors. As a demonstration of the analysis method these have 
been extracted using the MCMC (over the four main parameters only) and all have been 
recovered with their known injection parameter values. The extracted parameters for the 
two strongest signals are shown in Figure Ei From these injections it can be seen that the 
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Fig. 2. — The extracted posteriors for the pulsar parameters from the two strongest hardware 
injections for each detector and for a joint detector analysis. In each plot the injected 
parameter value is marked by a dashed vertical line. The left pulsar injection was at a 
frequency of 108.9 Hz and the right pulsar was at 194.3 Hz. 



parameters can be extracted to a high accuracy and consistently between detectors. The 
offsets in the extracted values from the injected values are within a few percent. These are 
well within the expected uncertainties of the detector calibration, which are approximately 
10%, 10% and 13% in amplitude, and 4?3 (0.08 rads), 3?4 (0.06 rads) and 2?3 (0.04 rads) in 
phase for HI, H2 and LI respectively. 



3. Evaluation of pulsar parameter errors 



For the majority of pulsars we have parameter correlation matrices that have been 
produced during the fit to the radio data, and (as discussed in §2.2p we can use these to 
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search over the uncertainties on the phase parameters. For some pulsars no correlation matrix 
was produced with the radio observations, so for these we instead construct conservative 
correlation matrices. These assume no correlations between any parameters, except in the 
case of binary systems for which we assume a correlation of 1 between the angle and time of 
periastron. This gives a slightly conservative over-estimate of the parameter errors, but is 
still a useful approximation for our purposes. From these correlation matrices and the given 
parameter standard deviations we produce a covariance matrix for each pulsar. 

Using these covariance matrices we can also assess the potential phase mismatch that 
might occur between the true signal and the best fit signal used in the heterodyne due to 
errors in the pulsar phase parameters (as discussed in §2.21 ) We take the mismatch (i.e. the 
loss in signal power caused by the heterodyne phase being offset, over time, from the true 
signal phase) to be 



where is the phase given a vector of phase parameters b at time t, and T is the total 
observation time. To get offsets in the phase parameters, 6b, we draw random values from 
the multivariate Gaussian defined by the covariance matrix. For each pulsar we drew 100 000 
random points and calculated the mean and maximum mismatch over the period of the S5 
run. The mean mismatch is a good indicator of whether the given best fit parameter values 
are adequate for our search, or whether the potential offset is large and a search including the 
phase parameters is entirely necessary. The maximum mismatch represents a potential worst 
case in which the true signal parameters are are at a point in parameter space approximately 
4.5(j from the best fit values (the maximum mismatch will obviously increase if one increases 
the number of randomly drawn values). Of our 113 non-glitching pulsars there are only 
16 pulsars with mean mismatches greater than 1%, three of which are greater than 10%: 
J0218+4232 at ~ 43%, J0024-7204H at ~ 15% and J1913+1011 at ~ 37%. Given the 
exclusion criterion for the S3 and S4 analyses (see §2.2p these three pulsars would have been 
removed from this analysis. There are 16 pulsars with maximum mismatches greater than 
10% with three of these being at almost 100% — i.e. if the signal parameters truly were offset 
from their best fit values by this much, and these parameters were not searched over, then 
the signal would be completely missed. This suggests that for the majority of pulsars the 
search over the four main parameters of h^, 0o, cost and is all that is necessary. But for 
a few, and in particular the three with large mean mismatches (J0218-I-4232, J0024— 7204H, 
and J1913+1011), the search over the extra parameters is needed to ensure not losing signal 





power. 
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4. Analysis 

We have used the MCMC search over all phase parameters, where they have given errors, 
for all but three pulsars (see below). For the majority of pulsars it is unnecessary (though 
harmless) to include these extra parameters in the search as their priors are so narrow, but 
it does provide an extra demonstration of the flexibility of the method. To double check, we 
also produced results for each pulsar using the four dimensional grid of the earlier analyses 
and found that, for pulsars with negligible mismatch, the results were consistent to within a 
few percent. 

For our full analysis we produced three independent MCMC chains with burn-in periods 
of 100 000 iterations, followed by 100 000 iterations to sample the posterior. The three chains 
were used to assess convergence using the tests discussed above in §2.11 All chains were seen 
to converge and were therefore combined to give a total chain length of 300 000 iterations 
from which the posteriors were generated. 

For the three pulsars that glitched during S5 (the Crab, J1952-I-3252 and J0539— 6910) 
the MCMC was not used to search over the position and frequency parameters as above, 
as the uncertainties on these parameters gave negligable potential mismatch. Our analysis 
of these pulsars did however include extra parameters in the MCMC to take into account 
potential uncertainties in the model caused by the ghtches. We analysed the data coherently 
over the full run as well as in stretches separated by the glitches that are treated separately. 
We also included a model that allowed for a fixed but unknown phase jump A(j) at the time 
of each glitch, keeping other physical parameters fixed across the glitch. Results for all of 
these cases are given in §5.11 

The maximum detector calibration uncertainties are approximately 10%, 10% and 13% 
in amplitude, and 4?3 (0.08 rads), 3?4 (0.06 rads) and 2?3 (0.04 rads) in phase for HI, H2 
and LI respectively. 



5. Results 

No evidence of a gravitational wave signal was seen for any of the pulsars. In light of 
this, we present joint 95% upper limits on ho for each pulsar (see Tabled] for the results for 
the non-glitching pulsars). We also interpret these as limits on the pulsar ellipticity, given 
by 

£ = 0.237 /i_24rkpc^^"'/38, (7) 
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where h-24 is the ho upper hmit in units of 1 x 10~ , I-ss = 1 and r^pc pulsar distance in 
kpc. For the majority of pulsars this distance i s taken as the estimate value from Australia 
Telescope National Facility Pulsar Cataloguj^ ( Manchester et al. 2005 ). but for others more 
up-to-date distances are kn own. For pulsars in T erzan 5 (those with the name Jl 748— 2446) a 
distance of 5.5 kpc is used (lOrtolani et al.ll2007l ). For J1939-I-2134 the best estim ate distance 
is a h ighly uncertain value of ~ 8.3 ± 5 kpc based on parallax measurements ( iKaspi et al. 



19941 ). but it i s thought to be a large overestimate, so instead we use a value of 3.55 kpc de- 
rived from the lCordes fc Lazid (120021 ) NE2001 galactic electron density model. The observed 
spin-down rate for globular cluster pulsars is contaminated by the accelerations within the 
cluster, which can lead to some seeming to spin-up as can be seen in column four of Table [1] 
For all pulsars in globular clusters, except J1824— 2452A, we instead calculate a conservative 
spin-down limit by assuming all pulsars have a characteristic age r = i>/2i> of 10^ years. Note 
that using the characteristic age gives a spin-down limit that is independent of frequency and 
only depends on r and r. Globular cluster pulsar J1824— 2452A has a large spin-down that 
is well above what could be masked by cluster accelerations. Therefore, for this pulsar we 
use its true spin-down for our limit calculation. For some nearby pulsars there is a small, but 
measurable, Shklovskii effect which will contaminate the observed spin-down. For these if a 
value of the intrinsic spin-down is known then this is used when calculating the spin-down 
limit. 

The results are plotted in histogram form in Figure [31 Also shown for comparison are 
the results of the previous search using combined data from the S3 and S4 science runs. 
The median upper limit on for this search, at 7.2x 10~^^, is about an order of magnitude 
better than in the previous analysis (6.3 x 10~^^). A large part of this increased sensitivity 
(about a factor of 4 to 5) is due the the longer observation time. The median ellipticity is 
e = 1.1x10"^, an improvement from 9.1x10"^, and the median ratio to the spin-down limit 
is 108, improved from 870. If one excludes the conservatively estimated spin-down limits for 
the globular cluster pulsars, the median ratio is 73. In Figure S] the upper limits on are 
also plotted overlaid onto an estimate of the search sensitivitj]^ . 
in these results due to the calibrations are given in §2.41 



The expected uncertainties 



The smallest upper limit on ho for any pulsar is 2.3x10 for J1603— 7202, which has 



http : //www. atnf . csiro . au/research/pulsar/psrcat/ 



^^The upper and lower estimated sensitivity limits for the shaded band in Figure |3] come from the values 
of the 95% ho upper limits that bound 95% of total values fro m simulations on white noise for randomly 
distributed pulsars. The band can be estimated from Figure 1 of lDupuis fc WoanI (|2005l ). which gives limits 
of (7 to 20) X ^y Sn/T, where S'„ is the single sided power spectral density and T is the total observation time 
in seconds. 
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a gravitational wave frequency of 135 Hz and is in the most sensitive part of the detectors' 
bands. The lowest ellipticity upper limit is 7.0 x 10"*^ for J2124— 3358, which has a gravita- 
tional wave frequency of 406 Hz and a best estimate distance of 0.2 kpc. Of the millisecond 
recycled pulsars this is also the closest to its spin-down limit, at a value of 9.4 times greater 
than this limit. Of all pulsars which did not glitch during S5, the young pulsar J1913+1011 
is the closest to its spin-down limit, at only 3.9 times greater than it. 



Table 1. Information on the non-glitching pulsars in our search, including the start and 
end times of the radio observations used in producing the parameter fits for our search, and 

upper limit results. 



Pulsar start - end (MJD) 


V (Hz) 


V (Hzs-i 


) 

J 


distance (kpc) 


spin-down limit 


joint hf°^" 


cllipticity 




t) / u 


10024— 72n4C^P 


48383 - 


54261 


173.71 


1.5 X 10" 


15 


4.9 


6 


.55 




10" 


-28 


5.88 


X 


10' 


-25 


2, 


.26 


X 


10" 


-5 


898 


10024— 72n4n'=P 


48465 - 


54261 


186.65 


1.2 X 10" 


16 


4.9 


6 


.55 


X 


10" 


-28 


4.45 


X 


10" 


-26 


1 


.48 


X 


10" 


6 


68 


10024— 72n4E'^^P 


48465 ~ 


54261 


282.78 


-7.9 X 10 


-15 


4.9 


6 


.55 




10" 


-28 


9.97 


X 


10 


-26 


1 


.44 


X 


10" 


6 


152 


70024- 72r)4F^P 


48465 - 


54261 


381.16 


-9.4 X 10 


-15 


4.9 


6 


.55 


X 


10" 


-28 


8.76 


X 


10" 


-26 


6, 


.98 


X 


10" 


-7 


134 


10024— 72040'=? 


48600 - 


54261 


247.50 


2.6 X 10" 


15 


4.9 


6 


.55 


X 


10" 


-28 


1.00 


X 


10' 


-25 


1 


.90 


X 


10" 


6 


153 


Tnn24— 72n4H'^'=P 


48518 - 


54261 


311.49 


1.8 X 10" 


16 


4.9 


6 


.55 


X 


10" 


-28 


6.44 


X 


10 


-26 


7, 


.69 


X 


10" 


-7 


98 


70024— 72047'^ '^P 


50684 - 


54261 


286.94 


3.8 X 10" 


15 


4.9 


6 


.55 


X 


10" 


-28 


5.19 


X 


IQ- 


-26 


7, 


.30 


X 


10" 


-7 


79 


70024— 7204 T'^'^p 


48383 - 


54261 


476.05 


2.2 X 10" 


15 


4.9 


6 


.55 


X 


10" 


-28 


1.04 


X 


10 


-25 


5, 


.34 


X 


10" 


-7 


159 


70024— 7204Tj'=p 


50687 - 


54261 


230.09 


6.5 X 10" 


15 


4.9 


6 


.55 


X 


10" 


-28 


5.82 


X 


10' 


-26 


1 


.27 


X 


10" 


6 


89 


70024— 72041X4^=? 


4849^ - 


54261 


271.99 


2.8 X 10" 


15 


4.9 


R 

U . 




X 


1 0" 


-28 


6.14 


X 


1 0" 


-26 


9 


61 


X 


1 0" 


-7 


94 


70024— 7204N'=P 


48516 - 


54261 


327.44 


2.4 X 10" 


15 


4.9 


6 


.55 


X 


10" 


-28 


8.35 


X 


10 


-26 


9, 


.02 


X 


10" 


-7 


128 


70024— 72040*^^? 


506Q0 - 


54261 


247.94 


-2.1 X 10 


-15 


4.9 


6 


.55 


X 


10" 


-28 


5.74 


X 


10' 


-26 


1 


.08 


X 


10" 


6 


88 


70024— 7204T! 


50743 - 


54261 


287.32 


-1.2 X 10 


-14 


4.9 


6 


.55 


X 


10" 


-28 


5.53 


X 


10' 


-26 


7, 


.76 


X 


10" 


-7 


84 


70024— 7204S'^'=P 


50687 - 


54241 


353.31 


1.5 X 10" 


14 


4.9 


6 


.55 


X 


10" 


-28 


6.82 


X 


10' 


-26 


6, 


.33 


X 


10" 


-7 


104 


TOnO/l 790/1 Tt>cp 


0U0o4 — 




131.78 


-5.1 X 10 


-15 


4.9 


D 


.00 


X 


1 n~ 

lU 


-28 




X 


1 D' 


-26 


9 


.Zo 


X 


1 n- 


6 


51 


J0024-72G4U^'=P 


48516 - 


54261 


230.26 


-5.0 X 10 


-15 


4.9 


6 


.55 


X 


10" 


-28 


5.63 


X 


10 


-26 


1, 


.23 


X 


10" 


6 


86 


J0024-7204Y''=P 


51504 - 


54261 


455.24 


7.3 X 10- 


15 


4.9 


6 


.55 


X 


10" 


-28 


9.42 


X 


10- 


-26 


5, 


.26 


X 


10" 


-7 


144 


J0218+4232^J 


49092 - 


54520 


430.46 


-1.4 X 10 


-14t 


5.8 


7, 


.91 


X 


10" 


-28 


1.47 


X 


10 


-25 


1, 


.10 


X 


10" 


6 


186 


J0407+1607^J 


52719 - 


54512 


38.91 


-1.2 X 10 


-16 


4.1 


3 


.48 


X 


10" 


-28 


6.18 


X 


10- 


-26 


3, 


.93 


X 


10" 


5 


178 


J0437-4715^P 


53683 - 


54388 


173.69 


-4.7 X 10 


-16t 


0.1 


8, 


.82 


X 


10" 


-27 


5.73 


X 


10- 


-25 


6, 


.74 


X 


10" 


-7 


65 


J0613-0200^^P 


53406 - 


54520 


326.60 


-9.8 X 10" 


-16t 


0.5 


2, 


.91 


X 


10" 


-27 


1.11 


X 


10 


-25 


1, 


.18 


X 


10" 


-7 


38 


J0621+1002^^ 


52571 - 


54516 


34.66 


-5.5 X 10 


-17t 


1.9 


5, 


.40 


X 


10" 


-28 


1.53 


X 


10- 


-25 


5, 


.65 


X 


10" 


5 


284 


JG711-6830P 


53687 - 


54388 


182.12 


-2.7 X 10 


-16t 


1.0 


9, 


.52 


X 


10" 


-28 


5.00 


X 


10 


-26 


3, 


.71 


X 


10" 


-7 


53 


J0737-3039A'^J 


53595 - 


54515 


44.05 


-3.4 X 10 


-15t 


1.1 


6 


.17 


X 


10" 


-27 


7.87 


X 


10- 


-26 


1, 


.10 


X 


10" 


5 


13 


J0751+1807^J 


53405 - 


54529 


287.46 


-6.3 X 10 


-16t 


0.6 


1, 


.92 


X 


10" 


-27 


1.64 


X 


10- 


-25 


2, 


.91 


X 


10" 


-7 


85 


J1012+5307^J 


53403 ~ 


54523 


190.27 


-4.7 X 10 


-16t 


0.5 


2, 


.43 


X 


10" 


-27 


6.94 


X 


10 


-26 


2, 


.36 


X 


10" 


-7 


29 


J1022+1001^^P 


53403 - 


54521 


60.78 


-1.6 X 10 


-16 


0.4 


3 


.27 


X 


10" 


-27 


4.44 


X 


10- 


-26 


1, 


.14 


X 


10" 


6 


14 


JlG24-0719jP 


53403 - 


54501 


193.72 


-6.9 X 10 


-16 


0.5 


2 


.88 


X 


10" 


-27 


5.01 


X 


10- 


-26 


1, 


.67 


X 


10" 


-7 


17 


J1045-4509*'P 


53688 - 


54386 


133.79 


-2.0 X 10 


-16t 


3.2 


3 


.00 


X 


10" 


-28 


4.37 


X 


10- 


-26 


1, 


.87 


X 


10" 


6 


145 
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Pulsar start - end (MJD) v (Hz) t'(Hzs ^) distance (kpc) spin-down limit joint /iq^'^" ellipticity h^'^" /h^ 



J1455- 


-3330^^ 


52688 - 


54524 


125.20 


-2.5 X 10 


-16t 


0, 


,7 


1 


.53 


X 


10" 


-27 


5, 


.15 


X 


10" 


-26 


5, 


.75 


X 


10" 


-7 


34 


J1600- 


-3053^P 


53688 - 


54386 


277.94 


-6.5 X 10 


-16t 


2, 


,7 


4, 


.62 


X 


10' 


-28 


5, 


.57 


X 


10" 


-26 


4, 


.55 


X 


10" 


-7 


121 


J1603- 


-7202^P 


53688 - 


54385 


67.38 


-5.9 X 10 


-17t 


1 


.6 


4, 


.62 


X 


10" 


-28 


2, 


.32 


X 


10" 


-26 


1 


.98 


X 


10" 


-6 


50 


J1623- 


-2631'"'=^ 


53403 - 


54517 


90.29 


-5.5 X 10 


-15 


2, 


,2 


1 


,46 


X 


10" 


-27 


5, 


.81 


X 


10" 


-26 


3 


.71 


X 


10" 


-6 


40 


J1640+2224^^ 


53410 - 


54506 


316.12 


-1.6 X 10 


-16t 


1 


,2 


4, 


.86 


X 


10" 


-28 


6 


.65 


X 


10" 


-26 


1 


.87 


X 


10" 


-7 


137 


J1643- 


~1224^JP 


52570 - 


54517 


216.37 


-6.8 X 10 


-16t 


4, 


.9 


2, 


.94 


X 


10" 


-28 


4, 


.35 


X 


10" 


-26 


1, 


.07 


X 


10" 


-6 


148 


J1701- 


-3006A'^'=P 


53590 - 


54391 


190.78 


4.8 X 10" 


15 


6 


.9 


4, 


.65 


X 


10" 


-28 


5, 


.82 


X 


10" 


-26 


2 


.61 


X 


10" 


-6 


125 


J1701- 


-3006B'^'=P 


53650 - 


54391 


278.25 


2.7 X 10" 


14 


6 


.9 


4, 


.65 


X 


10" 


-28 


7, 


.63 


X 


10" 


-26 


1, 


.61 


X 


10" 


-6 


164 


J1701- 


_3QQg(-;bcp 


53590 - 


54396 


131.36 


1.1 X 10" 


15 


6 


.9 


4, 


.65 


X 


10" 


-28 


3 


.52 


X 


10" 


-26 


3 


.32 


X 


10" 


-6 


76 


J1713^ 


f0747''^P 


53406 - 


54509 


218.81 


-3.8 X 10 


-16t 


1, 


,1 


9, 


.54 


X 


10" 


-28 


4, 


.44 


X 


10" 


-26 


2, 


.45 


X 


10" 


-7 


47 


J1730 


-2304jP 


52571 - 


54519 


123.11 


-3.1 X 10 


-16 


0, 


.5 


2, 


.49 


X 


10" 


-27 


5, 


.93 


X 


10" 


-26 


4, 


.72 


X 


10" 


-7 


24 


J1732- 


-5049^P 


53725 - 


54386 


188.23 


-4.9 X 10 


-16 


1 


.8 


7, 


.18 


X 


10" 


-28 


5, 


.25 


X 


10" 


-26 


6 


.34 


X 


10" 


-7 


73 


J1744- 


_1134jp 


52604 - 


54519 


245.43 


-4.1 X 10 


-16t 


0, 


.5 


2, 


.18 


X 


10" 


-27 


1 


.10 


X 


10" 


-25 


2, 


.07 


X 


10" 


-7 


50 


J1748- 


-2446A'^'=Jff 


52320 - 


54453 


86.48 


2.2 X 10" 


16 


5, 


.5 


5, 


.83 


X 


10" 


-28 


3 


.89 


X 


10" 


-26 


6 


.77 


X 


10" 


-6 


67 


J1748- 


- 24460^^^5 


53403 - 


54516 


118.54 


8.5 X 10" 


15 


5, 


.5 


5, 


.83 


X 


10" 


-28 


5, 


.00 


X 


10" 


-26 


4, 


.63 


X 


10" 


-6 


86 


J1748- 


-2446D=9 


50851 - 


53820 


212.13 


-5.7 X 10 


-15 


5, 


.5 


5, 


.83 


X 


10" 


-28 


6, 


.78 


X 


10" 


-26 


1, 


.96 


X 


10" 


-6 


116 


J1748- 


-2446E'^'=ff 


53193 - 


53820 


455.00 


3.8 X 10" 


15 


5, 


.5 


5, 


.83 


X 


10" 


-28 


8, 


.95 


X 


10" 


-26 


5, 


.62 


X 


10" 


-7 


153 


J1748- 


-2446F'=9 


53193 - 


53820 


180.50 


— 1.3 X 10 


-16 


5 


.5 


5 


.83 


X 


10" 


-28 


8 


.37 


X 


10" 


-26 


3 


.34 


X 


10" 


-6 


143 


J1748- 


-2446G'=f 


51884 - 


53820 


46.14 


-8.4 X 10 


-16 


5, 


.5 


5, 


.83 


X 


10" 


-28 


5, 


.82 


X 


10" 


-26 


3 


.56 


X 


10" 


-5 


100 


J1748- 


-2446^9 


51884 - 


53820 


203.01 


3.4 X 10" 


15 


5, 


.5 


5, 


.83 


X 


10" 


-28 


7, 


.81 


X 


10" 


-26 


2, 


,46 


X 


10" 


-6 


134 


J1748- 


-2446l'"'=ff 


50851 - 


54195 


104.49 


7.3 X 10" 


16 


5, 


.5 


5, 


.83 


X 


10" 


-28 


3 


.54 


X 


10" 


-26 


4, 


,21 


X 


10" 


-6 


61 


J1748- 


■ 2446X^=9 


51884 - 


53820 


336.74 


1.1 X 10" 


14 


5, 


.5 


5, 


.83 


X 


10" 


-28 


6 


.67 


X 


10" 


-26 


7, 


,65 


X 


10" 


-7 


114 


J1748- 


-2446L'=f 


51884 - 


53820 


445.49 


3.4 X 10" 


15 


5, 


.5 


5, 


.83 


X 


10" 


-28 


1, 


.39 


X 


10" 


-25 


9, 


,09 


X 


10" 


-7 


238 


J1748- 


2446M*^'=9 


51884 - 


53820 


280.15 


-3.9 X 10 


-14 


5, 


.5 


5, 


.83 


X 


10" 


-28 


1, 


.01 


X 


10" 


-25 


1, 


,68 


X 


10" 


-6 


173 


J1748- 


-2446N'^'=9 


53193 - 


54195 


115.38 


-7.4 X 10 


-15 


5, 


.5 


5, 


.83 


X 


10" 


-28 


5, 


.85 


X 


10" 


-26 


5, 


,71 


X 


10" 


-6 


100 


J1748- 


-24460'''=s 


52500 - 


53957 


596.43 


2.5 X 10" 


14 


5, 


.5 


5, 


.83 


X 


10" 


-28 


2, 


.65 


X 


10" 


-25 


9 


,68 


X 


10" 


-7 


454 


J1748-- 




53193 - 


54557 


578.50 


-8.7 X 10 


-14 


5, 


.5 


5, 


.83 


X 


10" 


-28 


1, 


.56 


X 


10" 


-25 


6 


,08 


X 


10" 


-7 


267 


J1748- 


■2446Q'''=9 


53193 - 


54139 


355.62 


4.6 X 10" 


15 


5, 


.5 


5, 


.83 


X 


10- 


-28 


8, 


.80 


X 


10" 


-26 


9, 


,05 


X 


10" 


-7 


151 


J1748- 


-2446R'=9 


52500 - 


53820 


198.86 


-1.9 X 10 


-14 


5, 


.5 


5, 


.83 


X 


10" 


-28 


8, 


.23 


X 


10" 


-26 


2, 


,71 


X 


10" 


6 


141 
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Pulsar start - end (MJD) v (Rz) ;> (Hzs ^) distance (kpc) spin-down limit joint ft,Q^^° ellipticity /iq^'^Y'^o'^ 



J1748- 


-2446S'=f 


53193 




53820 


163.49 


-1.7 X 10 


-15 


5, 


.5 


5, 


.83 


X 


10" 


-28 


4, 


.46 


X 


10" 


26 


2, 


.17 


X 


10" 


6 


76 


J1748- 


-2446T'=f 


51884 




53819 


141.15 


-6.1 X 10 


-15 


5, 


.5 


5, 


.83 


X 


10~ 


-28 


5, 


.12 


X 


10" 


-26 


3 


.34 


X 


10" 


-6 


88 


J1748- 


-2446V'''=s 


53193 




53820 


482.51 


2.2 X 10^ 


14 


5. 


.5 


5, 


.83 


X 


10" 


-28 


1, 


.26 


X 


10" 


-25 


7, 


.04 


X 


10" 


-7 


216 


J1748- 


2446W'^'=f 


52500 




53820 


237.80 


-7.1 X 10 


-15 


5, 


.5 


5, 


.83 


X 


10" 


-28 


9 


.57 


X 


10" 


26 


2, 


.20 


X 


10" 


-6 


164 


J1748- 


-2446X'''=s 


51884 




54139 


333.44 


-6.5 X 10 


-15 


5, 


.5 


5, 


.83 


X 


10" 


-28 


8, 


.18 


X 


10" 


26 


9 


.57 


X 
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Fig. 3. — The solid histograms show the results of this analysis in terms of upper limits on 
ho, the ellipticity e and the ratio to the spin-down limit (excluding the glitching pulsars). 
he clear histogram s show the same set of values for the combined S3 and S4 analysis 



flAbbott et all 120071 ). 



5.1. Glitching pulsars 



For the three pulsars that glitched, we have chosen to take into account three different 
models related to the coherence of the gravitational wave signal and the electromagnetic 
signal over the glitch: i) there is coherence between them over the glitch (i.e. the glitch 
causes no discontinuity between the electromagnetic and gravitational wave phases); ii) there 
is decoherence (in terms of a phase jump) between them at the time of the glitch, but the 
phase discontinuity is included as an extra search parameter (i.e. a A0 parameter is added 
at the time of the glitch); and iii) the data stretches before, between and after the glitches 
are treated separately and analysed independently. 



5.1.1. Crab pulsar 

For this search we include the three different models described above relating to the 
observed glitch in the pulsar on 2006 August 3. Potentially all the four main gravitational 
wave signal parameters could be changed during the glitch if it is large enough to cause major 
disruption to the star, but this is not the case for the observed glitch. It had a fractional 
frequency change of order Ai/fu ~ 5xl0~^, which is unlikely to be energetic enough to cause 
changes in the gravitational wave amplitude near our current levels of sensitivity. To model 



Table 1 — Continued 



Pulsar start - end (MJD) ly (Rz) j> (Hzs ^) distance (kpc) spin-down limit joint /iq^'^° ellipticity /iq^^°//iq'^ 



J1824-2452H^'=f 


53629 - 


54202 
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-3.6 X 10 


-15 


4, 


.9 


6, 
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-28 


8, 


.27 
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10 
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2, 
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X 
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6 
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J1824-2452J^'=9 


53629 ~ 


54201 
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4.7 X 10- 


15 


4, 


.9 


6. 


,55 


X 
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1, 
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X 


10^ 
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2, 


.03 


X 
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-6 
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.2 


4, 
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.65 
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'^The pulsar is within a binary system. 

'^Thc pulsar is within a globular cluster. 

SThe pulsar was observed by the Green Bank Telescope. 

JThe pulsar was observed by the Jodrell Bank Observatory. 

PThe pulsar was observed by the Parkes Observatory. 
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■^sensitivity estimate 

★ joint 95% upper limit (S5) 

• joint 95% upper limit (S3/S4) 
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Fig. 4. — The gravitational wave amplitude upper limits are plotted over the estimated 
sensitivity of the search as defined by the grey band (see text). Also plotted are the limits 
from the S3/S4 search. 
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the s i gnal p hase we use the regularly updated Crab pulsar Monthly Ephemeris flLyne et al. 



19931 . |2009[ ). which is needed to take into account the phase variations caused by timing 



noise. 



As for the previous Crab Pulsar search (lAbbott et al.ll2008l ) we have infor mation on the 



orientation of the pulsar from the orientation of the pulsar wind nebula (PWN) (INg &: Romani 



20081 ). We use this information to set Gaussian priors on the ip and l parameters of 
ip = 125?155 ± 1?355 (the dependence wraps around at ±45°, so the actual value used 
is 35?155) and l = 62?17 ± 2?20. 

There is reason to believe that the PWN orientation reflects that of the central pulsar, 
but in case this is not an accurate description we present results using both a uniform prior 
over all parameters, and using the restricted prior ranges. The MCMC only searches over 
the four main parameters and does not include errors on the pulsar position, frequency 
or frequency derivatives as these are negligible. The results are summarised in Table [3l 
where for the ellipticity and spin-down limit calculations a distance of 2kpc was used. The 
orientation angle suggested by the restricted priors is slightly favourable in terms of the 
observable gravitational wave emission it would produce, and this allows us to set a better 
upper limit. 

In Figure [5] we plot the result from model i), with restricted priors, as an exclusion region 
on the moment of inertia-ellipticity plane. It can been seen that for all the allowed regions 
in moment of inertia we beat the spin-down limit. If one assumes the Crab pulsar to have 
a moment of inertia at the upper end of the allowed range our result beats the spin-down 
limit by over an order of magnitude. 



5.1.2. PSRJ0537-6910 



For pulsar J0537— 6910, which has so far been timed only in X-rays, we rely on data 
from the Rossi X-ray Timing Explorer (RXTE) satellite. This pulsar is young, has a high 
spin-down rate, and is a prolific glitcher, and the refore we require o b servat ions overlapping 
with our data to produce a coherent template. iMiddleditch et al.l ( l2006l ) have published 
observations covering from the beginning of S5 up to 2006 August 21, during which time 
the pulsar was seen to glitch three times. Further observations have been made which 
span the rest of the S5 run and show another three glitches during this time. The epochs 
and parameters for the seven ephemeris periods overlapping with our data run are given in 
Table [21 For the first epoch there was no data for LI, so the joint result only uses HI and 
H2 data. For the analyses using all the data (models i and ii) we have 474 days of HI data. 
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475 days of H2 data and 397 days of LI data. Due to the glitches we perform parameter 
estimation for the same three models given above. 

As with the Crab Pulsar there is also information on the orientation of J0537— 6910 from 



model fits to its pulsar wind nebula (INg &: Romanil 120081 ). These are used to set Gaussian 



priors on ip and l of ip = 131?0 ± 2?2 (equivalently 41?0) and l = 92?8 ± 0?9. We again 
quote results using uniform priors over all parameters and with these restricted priors. The 
distance used in the ellipticity and spin-down limits for J0537— 6910 is 49.4 kpc. The results 
are summarised in Table [3l Using the restricted priors we obtain a worse upper limit on 
ho than for uniform priors. This is because the nebula suggests that the star has its spin 
axis perpendicular to the line of sight, and therefore the gravitational radiation is linearly 
polarised and the numerical strain amplitude is lower than average for a given strain tensor 
amplitude /lo- 
in Figure [5] we again plot the result from model i), with restricted priors, on the moment 
of inertia-ellipticity plane. It can be seen that the spin-down limit is beaten if we assume a 
moment of inertia to be 2 x lO^^kgm^ or greater. 



5.1.3. PSRJ1952+3252 

PSR J1952+3252 is another young pulsar with a high spin-down rate (although a couple 
of orders of magnitude less than for the Crab pulsar and J0537— 6910) - it has a spin 
parameters of z/ = 25.30 Hz and v = — 3.73 x 10~^^ Hz s~^. Jodrell Bank observations of 
this pulsar were made over the whole of S5, but it was observed to glitch at some point 
between 2007 January 1 and January 12. For both the pre and post-glitch epochs we have 
coherent timing solutions and again perform analyses as above. The are no constraints on 
the orientation of this pulsar, so we do not use any restricted priors. The results for this 
pulsar are given in Table O We reach about a factor of two above the spin-down limit using 
a distance of 2.5 kpc. The result from model i) is also plotted on the moment of inertia- 
ellipticity plane in Figure O It can be seen that I^s just over over 4, which is above the 
expected maximum allowable value, would be needed to beat the spin-down limit. 



6. Conclusions 

In this paper we have searched for continuous gravitational waves from an unprece- 
dented number of pulsars with unprecedented sensitivity, using coincident electromagnetic 
observations of many millisecond and young pulsars. Our direct upper limits have beaten the 



Table 2. Ephemeris informatio n for PSR J0537— 6 9 10. T he first four values are taken from 

Middleditch et all tood ). 





start - end (MJD) 


(Hz) 


V (Hzs-i) 




i) (Hzs 




epoch (MJD) 


1. 


53551- 


-53687 


62.000663106 


-1.994517x10" 


10 


11.2x10 


-21 


53557.044381976239 


2. 


53711- 


-53859 


61.996292178 


-1.993782x10- 


10 


8.1x10- 


-21 


53812.224185839386 


3. 
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61.995120229 


-1.994544x10- 
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-21 


53881.096123033291 


4. 


53953- 


-53996 


61.993888026 


-1.995140x10" 


10 


9.6x10- 


-21 


53952.687378384607 


5. 


54003- 


-54088 


61.992869785 


-1.994516x10" 


10 


2.2 x 10" 


-21 


54013.061576594146 


6. 


54116- 


-54273 


61.990885506 


-1.994577x10- 


10 


8.1x10- 
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Fig. 5. — The results of the Crab pulsar, J0537— 6910 and J1952-f3252 analyses, and the 
spin-down limits, plotted on the moment of inertia-ellipticity plane. The results used are 
those from model i) and with restricted priors on the angular parameters for the Crab pulsar 
and J0537— 6910. Areas to the right of the diagonal lines are excluded. The shaded regions 
are those outside the theoretically predicted range of moments of inertia Jss =1-3. 
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Table 3. Results of the analysis for the Crab pulsar, J0537— 6910 and J1952+3252. 



ellipticity hf'^''/hf 
epoch uniform restricted^ uniform restrictedt uniform restricted^ 



Crab pulsar 



model i)^ 


2.6x10 


-25 


2.0x10- 


-25 


1.4x10-"^ 


1.1x10 


-4 


0.18 


0.14 


model ii)*^ 


2.4x10 


-25 


1.9x10- 


-25 


1.3x10-^^ 


9.9x10 


-5 


0.17 


0.13 


1. 


4.9x10 


-25 


3.9x10- 


-25 


2.6x10-4 


2.1x10 


-4 


0.34 


0.27 


2. 


2.4x10 


-25 


1.9x10- 


-25 


1.3x10-4 


1.0x10" 


-4 


0.15 


0.13 


J0537-6910 


model i)^ 


2.9x10 


-26 


3.9x10- 


-26 


8.9x10-^ 


1.2x10" 


-4 


1.0 


1.3 


model ii)*^ 


4.1x10^ 


-26 


4.6x10" 


-26 


1.2x10^4 


1.4x10" 


-4 


1.4 


1.5 


1. 


9.1x10- 


-25 


1.2x10- 


-24 


2.8x10^3 


3.7x10" 


-3 


31.2 


41.2 


2. 


1.2x10" 


-25 


1.5x10- 


-25 


3.6x10-4 


4.6x10" 


-4 


4.1 


5.2 


3. 


1.4x10 


-25 


1.3x10- 


25 


4.2x10-4 


3.8x10 


-4 


4.7 


4.3 


4. 


8.3x10 


-26 


1.1x10- 


-25 


2.5x10-4 


3.3x10 


-4 


2.8 


3.7 


5. 


1.4x10" 


-25 


1.3x10- 


-25 


4.1x10-4 


4.0x10" 


-4 


4.7 


4.6 


6. 


4.4x10" 


-26 


5.8x10- 


26 


1.3x10-4 


1.8x10" 


-4 


1.5 


2.0 


7. 


4.9x10" 


-26 


6.1x10- 


26 


1.5x10-4 


1.9x10" 


-4 


1.7 


2.1 


J1952+3252 


model i)^ 


2.5x10" 


-25 






2.3x10-4 






2.0 




model ii)*^ 


3.6x10 


-25 






3.3x10-4 






2.9 




1. 


4.7x10 


-25 






4.3x10-4 






3.8 




2. 


4.4x10" 


-25 






4.0x10-4 






3.5 





^Uses observationally restricted priors on the orientation angle and polarisation angle 
^Uses the full data set for a coherent analysis 

^Uses the full data set, but searches over an extra phase parameter at each glitch 
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indirect limits (including the spin-down limit) for the Crab Pulsar and are at the canonical 
spin-down limit for J0537— 6910. For several more pulsars our upper limits are within roughly 
one order of magnitude of the spin-down limits, and therefore we expect a comparable search 
of Advanced LIGO and Virgo data to beat the latter. 

For the Crab Pulsar we improve upon the previous upper limit, which already beat the 
spin-down limit and other indirect limits; and we also provide results for different scenarios 
in which the signals may not be fully coherent over the run. By assuming the observationally 
constraine d priors on the angul ar parameters are correct we can limit the power output of the 
pulsar (see lAbbott et al.l ( 120081 )) to be less than ~ 2% of the canonical spin-down luminosity. 



Pulsar J0537— 6910 is an interesting case due to its youth and glitchiness. Assuming 
that the gravitational wave signal phase is coherent with the X-ray pulses over glitches, our 
upper limit is within a few percent of the canonical spin-down limit for this pulsar. If we 
assume that the phase is incoherent over the glitches then using the full data set our upper 
limit is ~ 1.2 times the canonical spin-down limit. Even if we assume that none of the four 
main pulsar parameters remains coherent during a glitch (i.e. there is a major rearrangement 
in the star's structure) then the results from the final two epochs individually are still only 
approximately two times the spin-down limit. All of these direct upper limits are within 
the uncertainties of the spin-down limit, and (unlike the Crab) there is no firmly observed 
braking index and thus no substantially stricter indirect limit. 

Pulsar J1952+3252 gets close to its spin-down limit, but even taking into account mo- 
ment of inertia and distance uncertainties, probably would still just miss beating it. It does 
however make an exciting target for future runs. 

Our observational upper limits on the ellipticities of the Crab and J0537— 6910 are in the 
vicinity of 1 Elastic deformations of this magnitude have been pre dicted as sustainable 
(jOweru I2OO5I : iHaskell et al.l 120071 : iLinI 120071 : iKnippel fc Sedrakiaru 120091) not for normal neu- 



tron stars but only for exotic forms of crystalline quark matter (jXull2003l : iMannarelli et al. 



20071 ). Such ellipticities could also be sustained by internal toroidal magnetic fields of order 
10^^ G depending on the field configuration, equation of state, and superconductivity o f 
the star JCutleilbood : lAkgun fc WassermanI [20071 : iHaskell et allboosl : IColaiuda et aPboOsh . 
Therefore our observations have achieved the sensitivity to detect some of the more extreme 
possibilities, and have constrained the internal magnetic fields of the Crab and J0537— 6910 
to be less than of order 10^^ G. However, it is important to be clear that the upper limits 
reported here do not constrain the properties of crystalline quark matter, because gravita- 
tional wave observations constrain t he true ellipt icity rather than the maximum ellipticity. 
This distinction w as raade clear in lOwen ( 2005 ), but h as not been properly enforced in 
more recent work ( (ILiru 120071 : iKnippel fc Sedrakiarul2009l ). Also, while it is the case that, as 
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anticipated by iHaskell et al.l (120071 ). the upper bounds on ellipticity reported here for some 
pulsars, notably J2124-3358, have pushed into a regime where elastic strains might support 
such deformations, detection was not expected, as such large ellipticities would conflict with 
the spin- down limit. 



Recent work by iHorowitz fc Kadaul (120091 ) suggests that the breaking strain of neutron 
star crusts may be an order of magnitude higher (10"^) than the highest values (10^^) 
previously used in estimates of maximum sustainable ellipticities. Their simulations are 
strictly applicable only to the outer crust (i.e. no neutron drip), but since the reason for 
the high breaking strain is very generic — the extreme pressure simply crushes away defects 
that contribute to early fracture — it may apply to the inner crust (the major contributor to 
ellipticity) as well. In that case normal neutron stars could sustain ellipticities close to 10~^, 
which is beyond reach of the present search (for those stars where our results are at or near 
the spin-down limit) but would be accessible with data from advanced interferometers. If 
the high breaking strain holds for mixed phases in hybrid stars (which is not clear due to 
the increasing importance of the s trong interac tion at high densities), it would bring their 
maximum ellipticities estimated by lOwenI (120051 ) up an order of magnitude to about the 10~^ 
values achieved here. 

In 2009 July the LIGO 4 km detectors (featuring some upgraded systems and titled 
Enhanced LIGO) and the upgraded Virgo detector began their sixth and second science 
runs respectively (S6 and VSR2). These upgrades aim to provide significantly better strain 
sensitivities than S5/VSR1. For LIGO these improvements will primarily be at frequencies 
greater than 40 Hz, but Virgo will improve at lower frequencies too and should outperform 
LIGO below about 40 Hz. LIGO and Virgo will closely approach and could potentially beat 
the spin-down limits not only for the Crab and J0537-6910 but for six more known pul- 
sars: J1952+ 3252, J1913+1011 J0737-3039A, J0 437-4715H, and the recen tly discovered 
J1747-2809 JCamilo et al.ll2009l ) and J1813-1749 dOotthelf fc Halpernllioogj ). Below 40 Hz 
Virgo could beat the spin-down limits for three more pulsars: the Vela pulsar, J0205-I-6449 
and J1833— 1034. Most of these eleven pulsars are young and may be prone to large levels 
of timing noise and glitches. Therefore to reach the best sensitivities it is essential to have 
radio and X-ray observations of these sources made in coincidence with the S6 /VSR2 run. 
In particular the extension of the RXTE satellite is crucial to our ability to perform targeted 
searches of J0537— 6910 during this time of enhanced detector performance. 

On a time scale of a few years, the Advanced LIGO and Virgo interferometers are 



^''^This pulsar sits on a strong LIGO spectral noise line caused by the violin mode resonance of the sus- 
pension wires, so the required sensitivity may only be reachable with Virgo. 
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expected to achieve strain sensitivities more than 10 times better for the pulsars in the 
current band, and will extend that band downward to approximately 10 Hz. Searches at 
such sensitivities will beat the spin-down limits from dozens of known pulsars and also enter 
the range of ellipticities predicted for normal neutron stars, improving the prospects for 
direct gravitational wave observations from these objects. 
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